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ABSTRACT 
The heat transfer and pressure drop characteristics of water flow- 
' ing in a vertical annulus under counter flow conditions were investigat- 
ed, The heat exchanger used consisted of a straight tube in tube test 
section with water as the heat transfer medium in both tubes. Convec- 
tion coefficients of heat transfer between the inner tube and the fluid 
flowing in the annulus were determined from the measured inner tube wall 
and average fluid bulk temperature and the energy transferred to the 
water flowing in the annulus, 

The friction factor data indicate that the proper choice of dia- 
meter to be employed in evahuating the Reynolds number is the equi- 
valent diameter, de = (d,- de and show that the transition from 
laminar to ‘turbulent flow occurs at Re #% 2000 as it does for flow in- 
Side tubes, The results for the laminar region are in excellent agree- 
ment with the theoretical predictions. Data in the turbulent region 
correlate well with that in the literature for commercial pipes. 

The heat transfer data cover a range of Reynolds numbers from 200 
to 5000. Two different types of runs were made~ equal capacitance, 
(equal mass rate of flow in both tubes for this particular investiga- 
tion) and constant wall temperature. The data indicate that the bulk 
temperature measurements were not reliable. However, since the measured 
temperature differences were reliable, it was possible to correlate both 
laminar and turbulent flow data by one empirical equation. The bulk tem- 
perature distribution for ce constant wall temperature run‘ along the 
annulus length was found to have a large effect on the heat transfer 


coefficients. It is not linear and should be taken into account. The 





data were best correlated by evaluating fluid properties at an average 
film temperature rather than an average bulk temperature as recommend- 


ed in most of the literature, 
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I. Introduction. 

The problem of heat transfer in annular spaces has not been ade- 
quately solved and is complicated by the many forms in which the en- 
pirical correlation equations are presented, eS] ; [4] , [5] .* 
Although a number of equations relating convection heat transfer co- 
efficients to system parameters have been presented in recent years, 

[ 1 | ; [2] » these provide correlations only for selected ranges 
of the pertinent variables, The development of proper relations for 
predicting annular heat transfer is made difficult, not only by an 
almost complete lack of data in the laminar region [14] , but also by 
insufficient data covering Both laminar and turbulent flow in the same 
annulus. 

It has been the objective of this investigation to obtain (1) fric- 
tion data for water flowing in the annulus of a vertical double tube 
heat exchanger over a range of Reynolds S saver 180 to 8000 and (2) heat 
transfer data for the same annulus over a range of Reynolds number 200 
to 5000 which was the limit of the equipment. The purpose of the fri- 
ction data was to determine the proper characteristic length to be 
utilized in evaluating the Reynolds number. The heat transfer results 
were treated in a manner to bring out their comparison with existing 
empirical correlations, and to establish a correlation which would cover 


the transition Reynolds number range. 


(*) Numbers in brackets refer to the Bibliography on page 19. 





Il, Description of Equipment. 

Heat Exchanger. The heat exchanger proper consisted of a single 
vertical 1,125 inch 0. D. x 0.995 inch I. D. 99.9% copper tube inside 
a 1.372 inch O. D. x 1.242 inch I. D. 99.9% copper tube 142 inches 
long. The annular diametric clearance was 0.117 inches and was main- 
tained by a helical arrangement of 36 stainless steel ball bearings, 
0.093 inch diameter, imbedded in the inner tube wall and fastened by 
solder. Details of construction are shown in Figures 1 and essential 
dimensions in Table l] 

Certain aspects are worthy of comment. The annular space was made 
sufficiently small so that reasonable annulus water velocities might be 
ama at low mass rates of flow. Inlet and outlet header spaces were 
made sufficiently large to permit proper measurement of the bulk tempera- 
tures. Copper-constantan immersion thermocouples were used to measure 
the inlet and outlet temperatures of both streams. A switch was built 
which enabled both absolute temperatures and temperature differences to 
be read. A schematic wiring diagram of the switching arrangement is 
shown in Figure 2. The circuit was designed to eliminate the interference 
from the circulating circuits existing between all the thermocouples. 

On the outer annulus tube 36 copper-constantan thermocouples were anchor- 
ed by soldering in 1/32 inch holes drilled half way through the tube 
wall, These were placed in groups of three, equally spaced around a 
diameter, at 11 inch intervals down the length of the tube. The in- 

side wall temperature of the inner tube was measured by means of a 
"traveling" thermocouple. Several designs were used and Figures 3 and 


4% show the details of construction. A Victory Corporation 2000 ohm 





glass enclosed thermistor was used for one design of the "traveling" 
thermocouple. It was calibrated in water with all leads in place and 
the calibration curve is shown in Figure 15. 

Pressure taps were drilled in the outer tube at a distance of 21 
inches from each end giving a distance of 100 inches between taps. Both 
mercury and water manometers were used to measure pressure drop. For 
reading the small pressure differentials a Fisher cathetometer was used 
which read to the nearest 0.01 centimeter. 

Auxiliary Equipment. The overall test assembly is shown in Figure 
5. A large hot water storage tank provided a heat reservoir of uniform 
temperature. The water was not changed during the investigation, there- 
fore assuring a minimum of entrapped gasses and particles in the circuit. 
Fischer and Porter ee were used to measure the flow rates 
for both hot and cold systems. These also served to indicate stability 
of the flow. The 36 outer wall thermocouple readings were taken with a 
Minneapolis - Honeywell 48 position recording potentiometer. These read- 
ings were recorded to a precision of + 0.2°F. All other thermocouples 
were read with a portable precision hand balance Rubicon potentiometer 
which was accurate to the nearest + 0.1 °F. The thermistor resistance 
was read with a wheatstone bridge balancing arrangement and was accurate 
to the nearest + 0.1 °F. 

The hedt exchanger proper was wrapped with 2 inch width asbestos 


cloth and then covered with 1 1/2 inch thick pipe lagging. 





III. Experimental Procedure. 

Friction: To obtain isothermal friction data, pressure drop 
measurements were made of the fluid flowing in the annular space at 
various flow rates. The pressure taps were so located to give an entrance 
length of 180 equivalent diameters. Pressure differences of less than 
50 inches of water were determined with a water manometer using a catheto- 
meter to secure accurate readings. For pressure differences of 50 inches 
of water or more a mercury manometer was used. The water temperature 
for most of the runs was essentially room temperature. Five check runs 
were made at a temperature of 100°F. For low flow rates an electric 
timer and a weighing tank were used to determine the mass rate of flow. 
These measurements were used to check the calibration of the flow meters, 
At low flow rates the flow meter meaSurements varied as much 6% with 
the flow rate obtained by weighing, but the overall variation was less 
than + 2%. 

Heat Transfer: The procedure for the equal capacitance 1uns and the 
constant wall temperature runs was essentially the same. For equal capa- 
citance, the mass flow rates in the annulus and inner tube were kept the 
game Since the specific heat for water was constant over the temperature 
range employed. Flow rates for the annulus were varied to cover a 
Reynolds number range of 200 to 5000. 

For constant wall temperature runs the inner tube mass flow rate 
was held constant at as high a value as permitted by the apparatus, 
corresponding to a Reynolds number of about 11,000 for all runs. For 
the annulus, the mass flow rate was varied to correspond to a range of 


Reynolds number of 200 to 2800. 





The hot water temperature was controlled by rheostats on the 
heaters and the energy input measured by a voltmeter and ammeter. 

For the hot water circulating system sufficient cold make-up water was 
added to keep a constant water level in the reservoir. 

Data for the runs were taken after five to ten minutes of steady 
state operation of the system. Steady state conditions were denoted by 
the constancy of flow rates and temperatures and were usually reached 
within an hour after changes in operating conditions had been made. Re- 
peat data recordings were made of all temperatures except those of the 
"traveling" thermocouple. 


The inner tube wall temperatures were determined by moving the 


"traveling" thermocouple up or down the tube in 12 inch increments. 





IV. Experimental Results. 

The friction results are shown in Figure 6 and tabulated in Table 
2. The table also presents the temperatures, pressure differences and 
flow rates for the individual runs. 

The heat transfer results are shown in Figures 7, 8, 9, 10 and 17 
and tabulated in Table 3. The table also presents the temperatures, flow 
rates and all other data used to calculate the system parameters as well 
as the parameters themselves. Inner wall temperatures were not correct- 
ed for thermocouple error as calculations of the temperature difference 
introduced by the exposed wires, using Schneider's flat plate analogy 
Le] , showed this effect to be insignificant. 

The inner wall temperatures were corrected according to laminar 
or turbulent flow in the inner tube for the physical distance of the 
thermocouple from the inner wall. This distance was 0.02 inches and 
amounted to 4% of the tube radius. Figure 14 shows Martinelli’s [ 16) 
and Eckert and Drake's (is) temperature profiles for turbulent and 
laminar flow respectively. A value of 0.625 was used for turbulent 
flow calculations and a value of 0.15 for laminar flow calculations. 

A correction to account for the temperature drop across the inner tube 
wall was also made. The value of thermal conductivity, k, for pure 
copper of 220 BTU/hr. ft.°F was used. This value is probably slightly 
high as the tube is only 99.9% pure copper. 

The average bulk temperatures of the annular fluid for constant 
wall temperature runs were corrected to account for the temperature 
distribution with annulus length. The three thermocouple readings at 


each of the 12 outside annulus wall positions were averaged and a 





dimensionless temperature, ty,- tei /teg ~ te; , was computed for 
each position. Simpsons Rule was used to calculate an average bulk 
temperature for each run. The average inner annulus wall tempera- 
tures were determined by an arithmetic average of the end wall tempera- 
tures. 

For the equal capacitance runs both the inner and outer annulus wall 
temperature distributions were considered linear and the average wall 
and average bulk temperatures determined on that basis. 

The system parameters are based on a film value temperature, ty, 


defined as 


tr = n+ Coen /2 Cy 


The data were also correlated using bulk temperatures and an aver- 


age bulk temperature as a basis for calculating system parameters. 





V. Sources of Error, 


The temperature measurement technique offered the princizai chance 
for error and all numerical results were a direct consequence of these, 
To minimize these errors, all thermocouples for bulk temperature measure- 
ments were similar and all manufactured thermocouples came from the same 
spool of wire and were made by the same method. A circuit arrangement, 
Figure 2, was built to a circulating current effects. Measured 
differences rather than absolute temperature measurements were used 
whenever possible. All readings were made on the same potentiometer 


and referenced against the same standard cell. 





Figure A. 


As shown in Figure A, the sum of ATa,; + dt14 should equal the sum 
of At,; 4 at34 . These measured temperature differences were 
totaled and found to vary an average of 2 5% Tais indicates the circuit 
design eliminated the circulation errors. 

The sensitivity of the potentiometer was + .001 millivolt or 
approximately + 0.04°F. Conversion of millivolts to degrees Fahren- 


heit was done from Leeds and Northrup Conversion Tables for Thermocouples. 





Flow rate variation was a second potential source of error. The 
flow meter readings were checked during the friction test by weighing 
the mass of water for a given time. The flow meters are considered 
accurate to + 2%. At least five readings were taken for each run and 
they did not vary more than 1%. 

Further sources of error which were determined to be insignificant 
included (a) expansion or contraction of the annulus tubes with tempera- 
ture fluctuations, (b) variation of the temperature during the time re- 
quired for observation, (c) heat loss along the thermocouple wire and 
(d) the bi-metallic thermoelectric effect at the connector terminals. . 
The latter connectings were designed and located to insure that both 
eeenaciion’ for each thermocouple lead would remain at the same tempera- 
ture thus assuring no effect on the results from this extra junéeion. 

An uncertainty analysis of the above effects on the system para- 
meters indicates that the repeatability of experimental data should 


be within + 10%. 
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VI. Discussion of Results. 

Friction: The isothermal friction data shown in Figure 6 are plot- 
ted as friction factor versus Reynolds number, based on an equivalent 
diameter. The change from laminar to turbulent flow is seen to occur 
approximately at a Reynolds number of 2000 which corresponds to the 
critical value for transition of flow inside tubes [10]. As it is 
desirable for comparison purposes to have the critical value of Reynolds 
number which conforms with that of flow in tubes, it was decided the 
equivalent diameter was the preferred characteristic defining length. 

As shown by Lamb [7] , besides the usual relations for friction 


factor, 


- plejde , (2) 


and Reynolds number, 
Re = de G 4 (3) 


for flow of fluids in an annulus an additional geometry parameter, 


B (di/d.) = + atl ll ag : (4) 
1+ (di/da)"4 ate 2 ol 


is needed, The limits of D (s/c) being 


Bld/d2)=1 when d/d,=O 


10 





as for flow in tubes, and 


B Cd fd.) = Ls when di/d. = 1 


as for flow between parallel plates. The functional relationship 
between the parameters for laminar flow is theoretically found to be 
fs FE g (ab/ds). 
“ (5) 
For the annulus used in this investigation 
giving a theoretical prediction of the friction factor for laminar 
flow as 


Seca f (6) 
Re 


From F .gure 6 it is seen that the experimental values of the friction 
factor are within + 5% of the theoretical values over a Reynolds 
number range of 200 to 2000. 

In the turbulent flow region, (Re > 2000), the data were plotted 


¢ 


and correlated by the equation: 


f= 0.0366 (Rey $ (dr/da) F (7) 


These data plot 12% higher than the proposed equation of Knudsen and 
Katz (3s) of: 


—O, 25 


f= 0.070(Re ; 


(8) 


and 6.5% higher than the proposed Davis [9] relationship: 


ry 





~O 


> 
f= 0.055(Re) (e247! wAdy= 1)" | 


da/d : (9) 


The correct function of the Reynolds number to use appears to 
be ea. as the slope of the line in Figure 6 is -0.2. Equation 
(8) expresses the friction factor by taking the average of a number 
of investigators’ experimental results [8] and neglects to account for 
any function of diameter. Knudsen and Katz [8] predict a 35% varia- 
tion from equation (8), depending on the ratio of (da/d, ). Only 
for larger (da/d,) ratios, above 3, may equation (8) be used with 
small error. It is, therefore, recommended that the Davis equation be 
used, especially for small values of (da/d, ). It is noted that the 
turbulent data fall very close to the line for flow inside commercial 
pipes. 

Heat Transfer, Figure 7 is a plot of all the laminar flow data 
correlated by using the Chen, Hawkins and Solberg [11} equation based 
on the arithmetic mean bulk temperature of the fluid and the viscosity, 


a evaluated at the inner annulus wall temperature t,, 


AS Ss of 14 4 9 
Nu z= 1.02 (Re) (Pr) (Gr) (Abfu,) (de/L) — @a/dy) (10) 
The data fall an average of 30% below the empirical curve. 
For the turbulent equal capacitance runs, Monrad and Pelton’s [12] 
equation based on the average fluid bulk temperatures, 
9 133 3 
Nu = .020 (Re) (Pr) (di/d)3; Re > 10,000, (11) 
was used and a plot is presented in Figure 8. The turbulent data fall 
50% below the empirical curve. 


The lack of correlation between these data and the empirical equations 


is obvious. The main explanation lies in the fact that the measured bulk 


vz 





temperatures were not true temperatures, The mixing chambers did not 
perform properly and further, the inner pipe conducted heat into both 
annuli mixing chambers, which would also contribute to bulk temperature 
error. This is substantiated by comparing the hot and cold bulk tempera- 
ture readings at the ends of the annulus. The cold water inlet tempera- 
ture which should remain fairly constant rises and falls with a hot 
water temperature rise and fall. 

The hot water bulk temperature measurements are considered better 
than the cold water bulk temperature measurements because the thermo- 
couples could be placed nearer the center of the tube. The cold water 
bulk thermocouples could not be placed closer than the outside radius 
of the inner tube and, therefore, could be considerably influenced by 
conduction from the hot tube as well as not actually "seeing" the in- 
let or outlet cold water temperatures, The inner and outer annulus wall 
temperatures are considered good as well as the bulk temperature differ- 
ences, 

Another possible reason for the discrepancies with these correla- 
tions is due to the effect of natural convection. There are no corre- 
lations in the literature that cover the flow conditions for the particu- 
lar annulus investigated. In general, natural convection is expected to 
influence the heat transfer behavior in an amount of + 15% depending on 
its relation to the external forces C13] . For the Chen, Hawkins and 
Solberg correlations with the Grashof number raised only to the 0.05 power, 
practically all of the free convection effects are damped out. In the 
Monrad and Pelton correlation no free convection effects are included 


since, with Reynolds numbers in excess of 10,000, the inertia forces 


is 





completely override the buoyant forces. 

It was found that all of the data of this investigation, in- 
ciuding the transition Reynolds number range, could be correlated by 
the equation 


SPY : , 065" ; 
Nu = 0.093 (Re) (Pr) (er) ces” aya”! 


(12) 

a plot of which is presented in Figure 10. This equation correlated 
all data with an average variation of + 10% and a maximum variation of 
20%. The Sieder and Tate (gl viscosity correlation, (bj), is not 
included as it is recognized that the present data do not cover a suf- 
ficient variation in temperature gradient to test this relationship, 
and also previous investigations have not been sufficiently extensive 
to prove this effect [2] . The parameters of equation (12) are 
evaluated at a film temperature defined by equation (1). The use of a 
film temperature allows for variations of all physical properties and 
resulted in the best correlation by an empirical equation. 

For equal capacitance laminar runs, the inner wall temperature 
correction for thermocouple location was made using Martinelli'’s [16] 
turbulent temperature distribution rather than the laminar temperature 
distribution. This was done because it gave the best correlation of 
these data vith respect to constant wall temperature data and the turbu- 
lent equal capacitance data, 

Figures ll, 12, 13, and 16 are dimensionless plots of annulus wall 
temperatures versus annulus length, and show that for the flow conditions 
of this investigation the temperature distributions are definitely not 
linear, However, it is believed that the arithmetic mean of the fluid 
bulk temperature is a close approximation to the average fluid tempera- 


ture for equal capacitance runs. The arithmetic mean bulk temperature 


14 





for the constant wall temperature runs was found to vary as much as 
25% from the average bulk temperature as found by applying Simpson's 
Rule. It is recommended that the average bulk temperature be used to 
evaluate the heat transfer as well as the film temperature for evalua- 
ting fluid properties, The average bulk temperature for the constant 


wall temperature runs was found to be 


Coen = »673 Ate + Ever ) (13) 


and is recommended for use where knowledge of the bulk temperature dis- 
tribution is unknown, 

All of the laminar runs were also correlated using the Chen, Hawkins 
and Solberg equation (10) but basing the fluid properties on a film tem- 
perature and using the above average bulk temperature, equation (13). 
These data plot an average of 15% higher than predicted by eerie (10). 
A plot is shown in Figure 9. 


A final correlation was made for all runs based on the Colburm an- 


alogy [19] of 
St (Pa)?= (Re) 

(14) 
with the fluid properties evaluated at the film temperature, equation 
(1). These results are shown in Figure 17. While it was expected that 
the Colburn analogy would produce an acceptable curve for the fixed 
geometry of the system and counter flow fluid arrangement, it was indeed 
Startling to find the small scatter of the data and the straight line 
behavior exhibited. It can only be postulated from the correlation that 
for the Reynolds number range covered and the small temperature differences 


which were employed, natural convection effects for all of the runs must 


15 





have been of the same order of magnitude. Any attempt to introduce 

the Grashof number into the correlation merely scattered the data, thus 
indicating the free convection effects were negligible. The extrapola- 
ted conclusion can be only that for a certain range of diameter ratios 

and equivalent diameters, even for laminar flow conditions, the Grashof 


number is not a pertinent parameter. 


16 





VII. Conclusions and Recommendations. 

The friction data indicate that the proper Reynolds number to 
employ is deGf, , where de is the equivalent diameter of the annulus. 
The friction data in the laminar region check the theoretical relation 
of Lamb C7] and this should be employed for annuli in this region of flow. 

The heat transfer data, both laminar and turbulent, were correlated 
by equation (12) based on the average bulk temperature and the fluid pro- 
perties evaluated at a film temperature defined by equation (1). This 
equation is presented as being valid only for use in investigating this 
particular geometry annulus and for identical conditions of flow. It 
is also suggested that equation (12) is valid for both heating and cool- 
ing. 

There is to date no single empirical correlation which can account 
for all the possible varistionsof flow conditions and annuli geometry. 

If knowledge of the heat transfer within + 5% is needed, then a particu- 
lar correlation must be found in the literature that applies co the 
specific type of flow pattern and annuli geometry under investigation. 

It would be interesting to know exactly how natural convection 
effects the heat transfer in vertical annuli. This could be found by 
reversing the direction of flow in the annulus and tube and analyzing 
the effect. It has to be expected that for large Reynolds numbers and 
small Grashof numbers, the influence of free convection on the heat trars- 
fer can be neglected. At what point this occurs is Gane to be determined. 
The Grashof number varied as much as 200% over the range of Reynolds 
numbers investigated. The Grashof number enters into equation (10) to 


the 0.05 power only and as such varied less than 2% over the range of 


17 





Reynolds numbers investigated. However, if the Grashof number enters 
imto the correlation to the 0.065 power as in equation (12), it will 
produce a 12% variation over the range of Reynolds numbers investigated. 
This agrees with Eckert and Diaguila's L1€] prediction of a 10% varia- 
tion in the heat transfer due to natural convection effects. An account 
must be made as to how the buoyancy forces act in conjunction with the 
external forces. A variation of 10% in the Nusselt number would result 


from ignoring this fact. 
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Table I, Apparatus Dimensions 
Inner Tube: 1 inch I. D. x 20 foot Copper Tube 
Effective tube length - 142 inches - 11.82 feet 
Diameter (inside) - 0.995 inches - 0.0825 feet 
Diameter (outside) - 1.125 inches - 0.0948 feet 
Wall Thickness - 0.065 inches 
Cross-Sectional Area (inside) - 0.00498 square feet 
Circumference (outside) - 0.29437 ft. 
Effective Heat Transfer Surface (outside) - 3.48 square feet 
Jacket Tube: 
Diameter (inside)-1.242 inches - 0.1035 feet 
Diameter (outside) - 1.3’2 inches - 0.31943 feet 
Cross Sectional Area (inside)~ 0.3193 feet 
Circumference (inside)- 4.46 inches - 0.372 feet 
Annulus: 
Equivalent Diameter: 0.117 inches - 0.00976 feet 
Cross-sectional Area - 0.00151 feet 


Effective friction length - 100 inches - 8.33 feet 
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TABLE II 


Isothermal Friction Data 


ys rate € ai —~ : 3 U 5 Re 
o ¥ /0 

Wsec, °F = Om oh, © Resa 0 ft/sec 

0.0208 64 8.16 0.716 On221 0.108 187 
0.0258 65 Li.20 0.705 0.274 0.0921 236 
0.0361 63 16.55 0.725 0.383 0.0696 322 
0.0450 62 20.64 0.736 0.477 0.0572 395 
0.0548 64 25.16 0.716 0.582 0.0457 495 
0.0665 64 29.88 0.716 0.705. 0.0371 600 
0.0826 63 33.08 @ 0.725 0.876 0.0307 736 
0.0945 64 42.53 0.716 1.005 0.0259 853 
0.1122 62 Zi) 0.736 0.190 0.0229 985 
0.1293 62 58.7 0.736 be 2 0.0192 23D 
0.1367 64 62% 7 0.716 1.448 0.0184 1232 
0.1412 62 69.0 0.736 1.498 0.0189 1240 
0.1823 61 95.3 0.747 1.932 0.0157 1575 
0.3152 62 201 0.736 Soo 0.0113 2760 
0.371 62 283 0.736 3.94 0.0111 3260 
0.494 60 453 0. 739 Dla 0.0102 4200 
0.569 58 587 0.780 6.04 0.00995 4720 
0.674 57 800 0.791 7.15 0.00965 5500 
0.751 37 978 0.791 7.97 0.00951 6130 
0.0411 98 12395 0.460 0.438 0.0415 578 
0.0807 98 202 0.460 0.862 0.0209 1132 
0.145 99 46,2 0.460 1.550 0.0118 2030 
0.338 99 196 0.463 3.61 0.00927 4320 
0.409 99.5 276 0.463 4.36 0.00895 5690 
0.533 99.5 460 0.467 5.69 0.00883 7360 
0.576 99.5 574 0.467 6.16 0.00932 7960 
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